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> I ■ Abstract Atom chips are a promising candidate for a scalable architecture for 

l' quantum information processing provided a universal set of gates can be imple- 

, mented with high fidelity. The difficult part in achieving universality is the entan- 

gling two-qubit gate. We consider a Rydberg phase gate for two atoms trapped 
on a chip and employ optimal control theory to find the shortest gate that still 
yields a reasonable gate error. Our parameters correspond to a situation where the 
Rydberg blockade regime is not yet reached. We discuss the role of spontaneous 
. emission and the effect of noise from the chip surface on the atoms in the Rydberg 

state. 

. Keywords Optimal control • phase gate ■ Rydberg atoms • cavity quantum 



electrodynamics 



1 Introduction 

0^ ' Neutral trapped atoms allow for qubit encoding in metastable internal states; they 

CO I offer long decoherence times and can easily be controlled by optical and magnetic 

■ fields. This makes them a promising candidate for the realization of a quantum 

computer [1]. Since the atoms hardly interact with each other, the largest difficulty 
. encountered when building such a quantum computer is the implementation of an 

' entangling two-qubit operation. The clock speed is thus typically limited by the 

gate operation time for one two-qubit gate which, together with several single- 
qubit operations, provides universal quantum computing 12,. 
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A number of schemes to entangle qubits carried by neutral atoms have been 
proposed [3l[ll[5l[6ir7l[8l[9l [T0ll 1 1| . An estimate of the corresponding gate duration is 
obtained in terms of the inverse of the interaction strength promoting the entan- 
glement. Adiabatic time evolution is not essential; it can be avoided by employing 
optimal control theory |12II13IIT4II151I16] . In fact, optimal control theory is an ex- 
tremely versatile tool that allows for achieving a high-fidelity implementation of 
desired quantum tasks. Moreover, it can be used to determine what fundamentally 
limits fidelity and duration of a quantum operation ^17j . 

For a two-qubit gate with neutral atoms, the highest gate speeds can be ex- 
pected when the atoms are internally excited to expose them to long-range interac- 
tions, a prominent example being dipole-dipole forces between Rydberg atoms [6l 
I18II19|[20] . However, care needs to be taken if the atoms are resonantly excited into 
a state with long-range interactions. For the example of the Rydberg gate, this 
will happen, e.g., if a perfectly entangling gate is desired but the distance between 
the atoms is too large to reach the Rydberg blockade regime [21] . The forces that 
induce entanglement will then couple internal electronic and vibrational dynamics, 
and the motional state of the atoms will have changed after the gate operation, 
implying strong leakage out of the quantum register. In principle, the vibrational 
excitation energy can be absorbed by an external field determined for example by 
optimal control theory. However, this requires the gate duration to be long enough 
to resolve the motional energy levels [22 . For resonant excitation into states with 
long-range interaction, the minimum duration, or quantum speed limit, for an en- 
tangling gate is thus limited either by the inverse interaction strength or by the 
vibrational period of the trap, whichever one is longer j2 11122] . 

Trapping, controlling and entangling two atoms is just the very first step to- 
wards a neutral atom based quantum computer. A scalable architecture is required 
to assemble many qubits and carry out any meaningful computation. One possi- 
bility to achieve scalability is given by miniaturizing the tools for trapping and 
manipulating the atoms to fit on a micrometer-sized chip [23.,24,,25 . These mi- 
crotraps offer strong confinement (fast vibrational motion) that may help in fast 
entangling gates. However, the close proximity of the atoms to the surface of the 
chip introduces noise sources which might significantly compromise the gate per- 
formance. Here, we study whether the fastest possible entangling gates for neutral 
atoms known to date - optimal Rydberg phase gates at the quantum speed limit 
- can be implemented on an atom chip, despite the presence of noise. To this 
end, we present an optimal control analysis for a pulsed excitation scheme of two 
atoms, using trap parameters realistic for microtraps. The analysis includes the nu- 
clear motion of the atom pair and addresses the influence of spontaneous emission 
from electronically excited states. We also estimate the infiuence of electric fields 
(static and dynamic) that arise from surface contamination and from blackbody 
and thermal near field radiation. Attention is payed to the specific properties of 
the fragile, highly excited Rydberg states that show huge electric dipole moments. 
Our results suggest that a fast gate (on the scale of some 10 ns) is possible with 
errors on the level of 10~^. We identify options for reducing this error by another 
order of magnitude to reach the fault tolerance threshold. 
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2 Fast phase gates via optimal control 

2.1 Controlled phase gate via Rydberg interaction of neutral atoms 

We consider two rubidium atoms that are trapped with a typical distance of 4 /im 
between them. The atoms sit in the ground state of the trap which is assumed to 
be deep enough such that each well can be approximated by a harmonic potential. 
The qubit states are encoded in two hyperfine levels of the ground state, for exam- 
ple |0) = 15si/2,F = 2,Mf = 2), jl) = !5si/2,F = 1,Mf = 1). A controlled phase 
gate can be implemented by excitation to a Rydberg state where the two atoms 
are exposed to a long-range interaction [6]. In a recent experiment with two ru- 
bidium atoms held in optical tweezers [2^, the \r) = j58d3/2,^" = 3,Mp = 3) 
Rydberg state was accessed by a near-resonant two-photon transition via the 
\i) = |5p2^/2:^ = 2,Mp = 2) intermediate state. The Hamiltonian for a single 
trapped atom in the rotating-wave approximation reads 

H*''w= E u)oi®(tf + ^,Vf)) 

+ ® (ff + Virapif) + ^) + !'^)(^! ® (^Tf + Vlrapif) + | 

+ £M!) (^|0)(^| + |^)(0|) ® 1(f) + ^ (\i){r\ + \r}{^\^ ® 1(f) . (1) 

Here, f denotes the position operator of the atom, T the kinetic energy opera- 
tor, and il\{t) (A = R, B) the time-dependent Rabi frequencies of the red and 
blue lasers (wavelengths 795 nm and 474 nm for the chosen Rydberg state) . The 
Rabi frequencies are parametrized according to f2\{t) = i7j.o(tanhe>,(t) -|- l)/2 G 
[0, Slj.o], where e\{t) will be determined by optimal control. A denotes the detuning 
of the red laser with respect to the intermediate state. The two-photon detuning 
from the Rydberg level is given in terms of the Stark shift, S = {Hg g — q)/4A. 
The corresponding level scheme is depicted in Fig. [TJ The electronically excited 
states |i) and \r) decay via spontaneous emission, their lifetimes are 27.7 ns for \i) 
and 210/iS for |r) (see also Sec. 13. 3p . The total Hamiltonian for the two atoms is 
given by 

H = ^ 2 ® If, + U.i ® Ifi ® H2^^ + Hmf ^ (2) 

with the interaction term describing the long-range forces when both atoms are in 
the Rydberg state, 

»lnf = \rr){rr\® % (3) 

The parameter C3 = 3.433- 10® a.u. = 3230 MHz pm^ has been calculated in Ref. [26] 
for the Rydberg levels |r) and scales with the principal quantum number as [20] . 
For an interatomic distance of rp = 4/im which is realistic for an experimental 
implementation, this yields an estimated interaction energy of 50 MHz. In all other 
states, the interaction between atoms that are thus far apart is so weak that it 
can be neglected. The harmonic approximation for the trapping potentials permits 
us to integrate over the center-of-mass motion of the two atoms and to focus on 
a one-dimensional model along the internuclear axis. We take an effective trap 
frequency uj/2tv ^ 276 kHz (period 3.6 /is) for the interatomic distance coordinate 
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Fig. 1 One-atom level scheme for near-resonant two-photon excitation to a Rydberg state. 



f = ]ri — f2| (ground state size « 20 nm). This value is slightly larger than in the 
optical tweezer traps of Ref . [26; , and may also be achieved in magnetic microtraps 
implemented on an atom chip [25] • The Hamiltonian, Eq. ([2]), is represented on 
an equidistant Fourier grid extending for ±0.3 fim around rp = 4 /im. The time 
evolution generated by this Hamiltonian is obtained by solving the time-dependent 
Schrodinger with the Chebychev propagator [gj. Note that we do not attempt 
here to spatially resolve the Rydberg orbits of the excited electron in the trapping 
fields. The level structure in a magnetic quadrupole field, for example, has been 
addressed in Refs. [28ll29l. 



2.2 Optimal control for two-qubit gates 



Our goal is to implement the controlled phase gate on the qubit register basis 
{|00),|01),|10),|11)}, 



O = diag 



(e'M,l,l) 



0= diag 1,1,1 



(4) 



This can be achieved by finding suitable fields e = {efl(i), eB(i)} that drive the 
system evolution from time zero to time T such that U(r, 0; e) = O up to a global 
phase. The quality of the implementation is measured by the fidelity, F. A suitable 
choice for evaluating it is based on projecting the actual evolution U(r, Oje) onto 
the desired operation O I16J, 



F=l[Tr[6+U(r,0;£) 



(5) 



A'' denotes the number of basis states. A*' = 4 if O is a two-qubit gate. Note that 
the actual evolution can proceed in a Hilbert space which is much larger than 
the qubit register space. In our example, the dimension of the Hilbert space is 
4 X 4 X Ar where Ar is the number of grid points required to represent f. Then 
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U(T, 0; e) is obtained by a suitable projection onto the register space which in our 

example is spanned by |(/7q-' ) = |j) ® ® \ fo), = 0,1, where \(po) represents 
the ground state of the trap. 

Optimal control treats the fidelity F, F £ [0, 1], as a functional of the control 
fields e. Allowing for additional constraints such as minimization of the integrated 
pulse energy, leads to the total functional J, 



J = -F + 



£ g{s)dt, g{s) = [sit) - e,ef (t)]' , (6) 



which is to be minimized. Here a is a weight, S{t) enforces the pulse to switch 
on and off smoothly, and eref(i) is the control field at the previous optimization 
step. Variation of the functional J with respect to the evolving basis states and 
the control field yields a set of coupled optimization equations that are solved 
iteratively. We employ here a linear variant of Krotov's method in order to obtain 
a monotonically convergent optimization algorithm [16j . In particular, the update 
equation for the field is then given by 

Z\£(t)=£('=+^)(i)-£W(t) = ^^{^f |6+0+(t,r;eW)AO(t,0;£('=+i')j^f 

jj' 

(7) 

which matches at each time t the target states, 0\(Pq'' ), propagated backward in 
time under the old fields, e^''\ with the initial states, i.e., the basis states \ipQ^ ), 
propagated forward in time under the new fields e^'^^^-'. 

In our example of the Rydberg gate, the short lifetime of the intermediate 
state, implies presence of a serious loss channel. The time evolution is then 
not unitary anymore. For dissipation due to spontaneous emission, the dynamics 
can be modeled by a Markovian Liouville-von Neumann equation for the den- 
sity matrix with Lindblad dissipators. Optimal control theory can be adapted to 
such a situation [30','31] but the numerical effort increases substantially. If dissi- 
pation is not needed to achieve the control objective but it rather represents a 
nuisance, a simpler approach is given by suppressing population in the dissipative 
channel. This is achieved by adding an additional constraint in the objective func- 
tional maximizing the expectation value of the projector onto the allowed subspace 
^aiiow = (1 ~ ® (1 — |*){*|) ® Ifi the subspace of the total Hilbert space 

excluding the unstable state [S^ , 

J=-F+ [ g{e)dt- [ {Paiiou,)dt. (8) 

Provided the modified control objective is fulfilled, i.e., the dissipative channel is 
never populated, the dynamics can again be described by a Schrodinger equation. 
Including the state-dependent constraint does not modify the update equation for 
the field, Eq. (O, except for a source term in the equation for the backward propa- 
gation [32J . Such an inhomogeneous Schrodinger equation can be solved efficiently 
with a modified Chebychev propagator [33J . We test the success of this scheme by 
propagating the initial states under the optimal field, and adding the imaginary 
term —^h'yi\i){i\ to the Hamiltonian ([!]), describing exponential decay from the 
intermediate state 
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2.3 Limits for tlie duration of optimal gates 

In the following, we employ the optimization algorithm with the standard func- 
tional ([6} as well as with the functional (|8]) that minimizes the population in the 
unstable state. We are looking for the best compromise between high fidelity and 
fast operation. 

The minimum gate durations are identified by running the optimization for 
varying T. If there is no loss, i.e., if we neglect spontaneous emission, our system 
is completely controllable. In the limit of sufficiently large T, the fidelity F then 
approaches one and the gate error approaches zero arbitrarily close [22 . As T is 
decreased, the time spent in the Rydberg state in which the atoms interact might 
become too short to pick up a non-local phase of tt (time scale Trfirp/cs ^ 10 ns). 
This is the first and most obvious constraint that limits the gate time T. Secondly, 
when the atoms are in the Rydberg state (note that our interatomic distance is 
too large for the blockade regime to be reached, i.e., the atoms are resonantly ex- 
cited into |rr)), they get accelerated in the attractive potential. In principle, 
the optimization can find laser fields that absorb this vibrational energy. However, 
in order to find such laser fields, the algorithm needs to distinguish between the 
motional target state, the ground state of the trap, and higher excited states. The 
vibrational period of the trap (r; 3.6 /xs) may therefore also limit the gate oper- 
ation time. For shallow traps and a strong interaction, this second limit will be 
larger than the first one, and thus determine the minimum gate time T ^2/2^. If the 
interaction cs/r'^ acts over a time such that a phase vr is picked up, then a momen- 
tum of order 5p = Znh/rQ is transferred to the relative motion. The motional error 
scales with the ratio to the momentum width in the ground state, (hujm/ A)^^'^ , 
as 1 — exp[— (5p^/(2fia;m)] « Sp^ /{2huim), see also Ref. 0- This translates into the 
condition that the ground state be much smaller than the mean distance. 



which is well met for our parameters since the trap ground state has a size ~ 20 nm. 

In our calculations the trap is only used to determine the initial state. During 
the action of the laser fields, the trapping potential is set to zero which corresponds 
to the dipole trap being switched off in the experiment [26]. We thus prevent the 
optimization algorithm from finding solutions that, if necessary, absorb vibrational 
energy. Once vibrational dynamics start to play a role, this will show up as a 
decrease in fidelity. However, this loss of fidelity could be avoided by leaving the 
trap on during the gate, in other words, it is not caused by a fundamental limit. 
Note that the trap is only needed to define the desired motional state at the end 
of the gate for the electronic ground state which carries the qubit state. 

Finally, a third factor limiting the gate time T is due to spontaneous emission 
from the intermediate state. While it turns out that this limits the best possible 
fidelity in current implementations of the Rydberg gate [26], it is not a funda- 
mental limit since one could think of accessing the Rydberg states differently, 
avoiding near-resonant two-photon excitation via such a strongly decaying state. 
We therefore discuss gate times and errors for both cases, including and neglecting 
spontaneous emission for the sake of comparison. 




(9) 
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gate time T 


gate 


error e 




no loss 


with loss 


30ns 


0.294 


0.305 


40ns 


0.043 


0.221 


50ns 


0.003 


0.025 


60ns 


0.003 


0.021 


70ns 


0.004 


0.021 


80ns 


0.005 


0.021 



Table 1 Optimal gate errors as a function of gate time. The optimization is performed here 
using the constraint that the population of the decaying intermediate state \i) be minimized 
[control functional Eq.l[8ll]. The right column gives the gate error when the optimal pulse is 
applied and an exponential loss rate is assumed from the state \i) (lifetime 27.7ns). 



2.4 Minimizing spontaneous emission: adiabatic pulse sequence 

Table [T] reports the minimal gate error e = 1 — F for a number of gate times 
T. Optimization was performed for a non-local phase in the |11) state, utilizing 
the control functional ([8]), i.e., we avoid spontaneous decay by minimizing the 
population of the intermediate state, |i) = |5pj^/2,-F = 2,Mp = 2). Other noise 
sources are neglected yet (see estimates in Sec. [3] below). The calculation allows 
for maximum Rabi frequencies SIji q = Ob.q = 27r ■ 260 MHz, cf. Ref. [26], and 
assumes a detuning A = 2-k ■ 600 MHz from |i) (see Fig. [T|). A number of Nr = 200 
spatial grid points is taken for the relative nuclear motion. For gate durations 
larger than 50 ns, gate errors of a few percent are obtained. The error increases 
only slightly when exponential loss from the state \i) is taken into account (column 
"with loss" ) . The difference between the gate error with and without loss illustrates 
how well the condition of suppressing population in \i) is fulfilled. For comparison, 
if the additional constraint is not employed in the optimization, the gate errors 
increase by two orders of magnitude to 20% and even more when spontaneous 
emission is taken into account. The gate errors without exponential decay in \i) 
(central column of Table [1} illustrate when excitation of the motional degree of 
freedom start to play a role - for gate durations of 70 ns and larger, the error 
obtained without spontaneous emission is mainly due to excitation of higher trap 
states. Once spontaneous emission is included, this effect is, however, not dominant 
anymore. 

The population dynamics of each electronic state induced by the optimal fields 
is shown in Fig. [2l demonstrating that indeed the intermediate state is almost 
never populated. The overall picture suggested by the population dynamics is an 
adiabatic transfer to and from the Rydberg state in a double-STIRAP-like fashion. 
This is confirmed by inspecting the shapes of the optimal fields in Fig. [3l The blue 
field connecting \i) and \r) plays the role of the Stokes pulse and the red field 
connecting jO) and \i) that of the pump pulse in the first half of the time interval. 
These roles are reversed in the second half. The middle part of the interval, between 
20 ns and 50 ns is crucial for obtaining the non-local phase, as can be seen from the 
significant population of the \rr) state. The dynamics of the jOO) and the |01), |10) 
states are locked. This is due to the fact that we optimize for a non-local phase x 
in the |11) state. However, the |11) state does not couple to the laser fields. The 
non-local phase, which is given by x = "/"oo — 4>oi ^ (pw + 4>ii [22,, is thus achieved 
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10 20 30 40 50 60 70 

time ( ns ) 



Fig. 2 Minimizing the population in the unstable state population dynamics for the 
sixteen electronic states (|11) is not shown since the population is equal to one for all times). 
The radiative lifetime of the states \i) and \r) is described by a non-Hermitian Hamiltonian. 
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Fig. 3 Optimal shapes of the laser fields and corresponding spectra, employing the constraint 
to minimize population of \i). 



by 'coordinating' the time evolution of the remaining three two-qubit states which 
can be controlled by the laser fields. 
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gate time T 


gate error e 


20ns 


0.067 


25ns 


0.051 


30ns 


0.001 0.003 


35ns 


0.001 


40ns 


0.001 


45ns 


0.002 


50ns 


0.003 


60ns 


0.003 



Table 2 Optimal gate errors as a function of gate time, neglecting losses from both the 
intermediate state |i) and the Rydberg state |r). The shortest possible gate time to achieve a 
reasonable fidelity is then limited by the interaction in the excited state and the motion in the 
trap. The error increases to — > 0.003 (at 30 ns) when a lifetime of 20 fis for the Rydberg level 
is included, consistent with the estimates of Sec. 13.31 



2.5 Ignoring spontaneous emission: Rabi flop sequence 

If the Rydberg state could be accessed directly, gate error e and time T would be 
solely determined by the interaction strength in the Rydberg state and the trap 
frequency 22 . We examine this scenario by ignoring the instability of the interme- 
diate state \i) and applying the standard control functional, Eq. Table[2]hsts 
the minimal gate error e neglecting the spontaneous decay from the intermediate 
state for a number of gate times T. Optimization was performed for a non-local 
phase in the |00) state. Inspection of Table [2] reveals that the minimum gate op- 
eration time is about 30 ns with a gate error of the order of 10~^. This limit is 
imposed by the interaction strength in the Rydberg state \rr). The fact that the 
gate error increases as T is enlarged, indicates that excitation of motion in the trap 
starts to play a role. Including the trapping potential in the calculation for the 
qubit states during the gate allows motional excitation to be further reduced [21] . 

The population dynamics of each electronic state induced by the optimal flelds 
shown in Fig.|3]are clearly non-adiabatic. The intermediate state |i) is significantly 
populated at intermediate times since the optimization algorithm does not 'know' 
that \i) corresponds to a loss channel. The population of the Rydberg state in the 
middle of the time interval is larger in Fig. |3]than in Fig. [2l This explains why a 
smaller gate duration can be achieved. The optimal pulse shapes and spectra are 
shown in Fig. (5] The observed sequence of fast switches of both red and blue laser 
fields pumps population in a ladder-like fashion with a time constant equal to the 
duration of a vr-pulse at the maximally allowed Rabi frequency. Correspondingly, 
we observe stronger sidebands in the spectra of the optimal pulses (compare Figs. [3] 
andO. 

In summary, we find minimum gate durations of about 30 ns to 50 ns for a high- 
fidelity implementation of the Rydberg gate. Shorter gate operation times might 
be possible by employing higher lying Rydberg states which exhibit a stronger 
interaction. However, higher states will also be more sensitive to stray fields. While 
our model includes all relevant degrees of freedom and spontaneous emission from 
the excited states, noise sources that are inevitably present in any experimental 
setup need to be analysed in order to gain a full understanding of what fidelities 
can be achieved in an experiment based on optimized pulse shapes. The infiuence 
of noise sources is expected to be particularly detrimental for miniaturized setups 
such as an atom chip. 
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Fig. 4 Two-photon excitation scheme treating the intermediate state |i} as stable: population 
dynamics for the sixteen electronic states. (|11) is not shown since the population is equal to 
one for all times.) 
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Fig. 5 Optimal shapes of the laser fields and corresponding spectra, assuming the intermediate 
state |i) to be stable. 



3 Influence of noise 



The hyperfine ground-state levels |0), as defined in Sec. 12.11 provide a qubit 
robust with respect to noise, as has been discussed elsewhere [23ll24l[25] and demon- 
strated experimentally [34j. We therefore focus here on the specific sensitivity of 
the Rydberg state \r) to static and fluctuating fields typical for an atom chip en- 
vironment. Since the Rydberg levels are populated only over the duration of the 
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laser pulse (some 10 ns) and their angular momenta are relatively small (j < 5/2), 
we do not have to consider magnetic field noise that takes effect only on a time 
scale of 100ms or longer [25, . 



3.1 Surface impurities: DC Stark shifts 

Previous work has shown that the surfaces of atom chips get contaminated with 
adsorbed atoms that lead to randomly placed charges or dipoles [35 . As a typical 
scenario, consider an alkali atom on a metal that gives off its electron into the 
surface so that an upright dipole moment remains. At a distance z ~ 10 /im above 
it, this single impurity creates a tiny electric field of the order of 

E,^,{z)^^^=0.l5^, (10) 

where as is the Bohr radius. But if 10^ adatoms are distributed over an area of 
{Wfim)^, (corresponding to an average distance ~ 100 nm, i.e., much less than a 
monolayer), their electric fields roughly add up to 0(1 V/m). The corresponding 
Stark-Hamiltonian, 

= -A-Ei,„p(z), (11) 

has matrix elements of the order eagn^ Ei^p{z) ~ 6.6 MHz for n = 58, just one 
order of magnitude below the Rydberg interaction Hamiltonian (|3]). The experi- 
ments of Ref. [35 have actually detected electric fields up to 1 kV/m at a distance 
of about 10 /im. On an insulating surface, charges can be trapped and the fields 
would even be stronger. Note that the threshold for field ionization is of the order 
of 5kV/m for the n = 58 level [36]. 

The Rydberg level considered here, \r) = |nd3/2) {n = 58), is actually "pro- 
tected" from a linear Stark shift because its quantum defect [5^ « 1.34) lifts the 
degeneracy with the opposite parity states |np, f, ...), by an energy splitting of 
the order of An = l/n^a.u. ^ 1.125cm-i = 33.7GHz (see Fig. |6] (left)). For this 
reason, the Stark shift is quadratic in the field and inversely proportional to the 
detuning (< An) from the nearest level with opposite parity. The nearest levels 
57f and 59p give quadratic Stark shifts that partially cancel each other, leav- 
ing a polarizability for 58d of the order of 50 kHz(V/m)~^. The 56f state shows 
a much larger polarizability because the hydrogen manifold 56H = 56g, h, . . . is 
only about 0.02 cm" ^ (^ 600 MHz) away (Fig. |6] (right)) and its influence is not 
cancelled by another level. We estimate a quadratic Stark effect of the order of 
2...3MHz(V/m)-2. 

This means that fields above 1 V/m start to detune the so-called Forster reso- 
nance between the two-atom states |58d3/2 58d3y2) and |56f5y2 60px/2) whose en- 
ergy mismatch is only 7 MHz « 2.3 • 10~'*cm~^ (see Fig. |6] (right) and Refs. [371 
126] V One then loses the strong scaling of the Rydberg interaction that turns 
into the weaker van der Waals scaling. In addition, at the level of 10 V/m, the 
linear DC Stark shift of the hydrogen-like state manifold |56H) exceeds ^ 600 MHz. 
By an avoided crossing, the state |56f5/2) is then pushed down, and the Forster 
resonance is detuned. 
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Fig. 6 (left) Energy levels near the Rydberg state |r) = |58d3/2)- The strong Rydberg 
interaction arises from the near degeneracy of \rr) with the two-atom level 56f5/2 60pi/2 
("Forster resonance" 1371 . dashed lines). The fine structure is not resolved on the scale 
zi58 = 1/58^ a. u. = 33.7GHz. The thick gray lines (label nH) give the hydrogen-like state 
manifold beyond the nf states (angular momenta ( > 3). The left and right columns show even 
and odd orbital angular momentum states. 

(right) Energy levels of the two-atom system. The levels in the Forster resonance are connected 
by the dashed line (detuning 7 MHz 1261 ). At a distance rg ~ 4/jm, the two levels hybridize 
and split by « 50 MHz, as described by the two-atom Hamiltonian ((Sj. Note the difference in 
scale: the fine structure of the p and d levels is resolved (but not for the f levels). A quantum 
defect of order 10~^ splits off the level 56f from the hydrogen-like manifold 56H = 56g, h, . . .. 
The gray triangles illustrate the linear Stark splitting of 56H in a weak static electric field be- 
tween and IV/m (shifts up to ±60 MHz). At slightly higher fields, also the levels 57f 59pi/2 
and 56f60p3/2 would be pushed away (avoided crossings). The quadratic Stark shift of the 
other levels is not visible on this scale. 

The energy levels are calculated from quantum defect data collected in Refs. I38II39I . A precise 
localization of the Forster resonance (not attempted here) would require knowledge of the 
quantum defects at the 10"* level. 



3.2 Fluctuating fields; dephasing 

The one- and two-atom Rydberg states that we are considering do not feature 
a permanent electric dipole moment. This is due to the quantum defects that 
split them off the hydrogen-like manifolds of higher angular momentum states 
{I > 3). In weak fields, their Stark effect is therefore quadratic. We consider here 
the beating between a fluctuating field and a static impurity field Simp (2) at 
the level of IV/m, small enough not to perturb the dipole-dipole interactions, 
see previous section. The mixing of states with opposite parity is linear in the 
impurity field and translates into a dipole moment of the order of dji = q^j J5imp {z) 
where an is the quasi-static polarizability. For the level 58d3/2, the contributions 
to nearby levels above and below partially cancel, leading to the relatively small 
value Id^^l ~ Aea^- The level 56i^^2 is much more polarizable, as found above, and 
dfl ~ 200 ea^- Note that this is still much smaller than expected from the average 
size of a Rydberg atom (« a^n^). The fiuctuating phase shift due to the Rydberg 
polarizability an is then given by 

A<t){T) = - jdt ^Ei^p(z) ■ Eauct{z,t) , (12) 
ti 

where the time integral is evaluated over that part of the pulse that the atom 
effectively spends in the Rydberg level (see Figs. [2l 
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10^ 10* 10^" 10^^ 10'^ 
frequency [Hz] 

Fig. 7 Electric field noise spectrum (in (V/m)^/Hz) due to patch charge fluctuations (dot- 
dashed) and due to Johnson-Nyquist noise from a gold half-space (solid and dashed lines). Dis- 
tance fixed to z = 10 fim, temperature 300 K. Patch fluctuations are calculated from Ref. 1401 
and extrapolate ion trap data (in the MHz range), assuming a l/ui scaling. The Johnson- 
Nyquist noise is calculated along the lines of Ref. |41| . It includes blackbody radiation and 
free-space vacuum fluctuations and changes from a cj^/^ into a oj'^ power law near the typical 
Rydberg transition energy A^g (arrow). Solid (dashed) lines are for flelds parallel (perpendic- 
ular) to the surface, respectively. The right scale gives the expected dephasing rate for a static 
dipole moment 0(200 ea^) and a white spectrum. The arrow at l/27rT marks the Fourier- 
limited band width for quasi-static noise over an effective interaction time r = 10 ns. The 
contribution of quantum fluctuations was subtracted in the Planck spectrum (dotted curve). 



Let US assume an electric field noise spectrum with a scaling 

5M^,-) = ¥^f^f (13) 



that arises from "patch charge" fluctuations, as observed in experiments with 
miniaturized ion traps |4(JII42|HH] (exponent /3 « 0.7 ... 1). The l/z^ scaling actually 
only holds at heights z beyond a characteristic length scale for the patch size. 
This effect is taken into account in the plots shown, using the model of Ref. [30] . 
The thermal electric fields that originate from the motion of charges in the chip 
material (Johnson-Nyquist noise) have a lower noise spectrum compared to the 
patch charge model of Eq. (|13p j 4 111421133] , at least in the low-frequency range (up 
to a few MHz) relevant for dephasing. This is illustrated in Fig. [7] where the two 
noise spectra (at fixed distance 10 /im) are plotted vs frequency. 

The dephasing of the Rydberg state is now estimated by calculating the vari- 
ance of the phase shift (|12p . Provided the interaction time t = t2 — ti is large 
compared to the noise correlation time, the variance increases like 

dj^SsizQ, l/r) 



A<l}^{r)^T 2;i2(^/^q)4 cos(7r/3/2)r(2 /3) ' ^^"^^ 

where the noise spectrum is evaluated at oj = l/r, roughly the Fourier-limited band 
width of the pulse, and r{-) is Euler's gamma function. For room-temperature mi- 
croscopic ion traps [33|130|1351|36] , noise levels of (75 fim. 1 MHz) ~ 10"^^ {V/mf/Hz 
and /3 « 0.7. . .0.8 are typical. Note that a low-frequency cutoff must be applied 
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Fig. 8 Normalized cross-correlation of electric field noise due to patch charge fluctuations, 
evaluated at three different positions and normalized to the noise spectrum in one trap. If the 
correlation is weak, the relative phase of the Rydberg state in the two traps is randomized on 
the same time scale as the common phase. The correlations are calculated by generalizing the 
model of Ref. [40'. Any group of three curves corresponds to patch correlation lengths 0.5, 1, 
and 1.5 fim. 



for pure 1/f noise. Taking an effective time of r ^ 20 ns spent in the Rydberg 
level, one gets at z = 10 fim a dephasing rate of 



D^{56i) = A(j>'^{t)/r = 0(10^ s"^) (15) 



for the highly polarizable Rydberg state. This affects in particular the two-Rydberg 
state via its admixture of 56f 60p (Fig. [6]). For the phase gate, we get a sizable 
decoherence error e « 1 — (e"^'^*''^^) ^ ^D^t ~ 1% from the uncertainty in the two- 
atom phase 011. We recall that this number scales quadratically with the impurity 
field, assumed here to be 1 V/m. The other Rydberg levels involved show a much 
smaller dephasing, in particular D^{58d) = 0(400 s~^). This is insignificant over 
the effective pulse duration r, and therefore the one-atom phases (j)oi, 4>io are not 
compromised. 

One may ask the question to what extent the fields seen by the two atoms are 
correlated. The phase gate misses its target already, of course, if there are fiuctu- 
ations that are common to both atoms. A differential phase would mix the even 
and odd states !56f5/2 6OP1/2) ± |60pi/2 56f5/2) that are involved in the Forster 
resonance. This can be quantified from the cross-correlation spectrum of the patch 
fields, Se{A, B; lj), where A and B are the positions of the two traps. In Fig. [8l we 
plot the normalized cross-correlation for two atoms at the same distance z from 
the chip, but laterally separated by a length tq. The two atoms are subject to the 
same noise (no differential dephasing) if the normalized correlation is unity. One 
notes that the fields decorrelate on a scale given by the height above the surface. 
The correlation is still quite strong for the parameters we considered here, i.e. a 
separation ro = 4/im and height 10 fim. 
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distance [//m] 

Fig. 9 Electric field noise spectra (in (V/m)^/Hz) vs distance from the chip surface, evalu- 
ated at the level splitting Ar^s typical for Rydberg transitions, and 300 K. The patch charge 
spectrum is extrapolated to this frequency assuming a l/o; scaling; the patch correlation length 
is in the range f = 1 it 0.5 fim (shaded area). The solid and dashed curves "clean Au" give the 
noise spectrum of thermal near field radiation, the far field limit is shown by the dotted line. 
The skin depth <5 separates two power laws ^1 , and for z > A/2, the spectrum oscillates due 
to interference by reflection from the surface. In this limit, a perfect conductor approach yields 
good agreement. Within the wide plateau betwen 10 fim and 1 mm, a dipole perpendicular 
to the surface (±) is subject to near field noise about five times as strong as the free space 
(Planck) level. 

3.3 Thermal radiation: lifetime and AC Stark shift 

The radiative lifetime of an emitter is strongly modified in the vicinity of a macro- 
scopic body, as illustrated by the Purcell effect [37j . In addition, thermal radiation 
plays a significant role because the Bohr frequencies of Rydberg atoms are low, 
typically An ^ kgT. It is well known that this reduces the lifetime and coherence 
time of trapped particles, even to the level that magnetic dipole transitions become 
relevant [25]. We estimate here surface-enhanced radiative loss and discuss that 
despite the large polarizabilities, thermal radiation gives rise to AC Stark shifts 
that are overall small, leaving only the zero-temperature van der Waals shift. A 
similar question has been addressed in Refs. [48.i49 _. 

Spontaneous emission in free space is dominated by the decay into the lowest 
states because of the cubic scaling of the vacuum field spectrum. This leads to a 
lifetime of the Rydberg state of 210 /is at T = 50 and a decay probability, i.e., 
gate error, of the order of 10~^ given the Rydberg excitation time r « 20 ns. We 
therefore need to check that this error does not increase significantly in an atom 
chip environment. The contribution of blackbody radiation reduces the lifetime 
significantly (absorption and emission), to a total figure of roughly 90 /is for the 
Rydberg level 58d where the effective quantum number is n* ^ 58 — 1.34 [50j . 

The chip surface enhances the mode density of the electromagnetic field [481 
1511152] , This leads to a different scenario, however, for transitions with large or 
small Bohr frequencies. For the decay into low states, the resonant wavelengths 
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are typically in the visible and near UV, small compared to the atom-surface dis- 
tance. These field modes form a shallow interference pattern due to reflection at 
the surface, enhancing or suppressing the decay by roughly a factor of 2. Destruc- 
tive interference can be used to suppress certain decay channels, as suggested in 
Ref . [53] in an application of the Purcell effect [47_ . 

At smaller Bohr frequencies, thermally stimulated emission and absorption are 
enhanced much more strongly in the near field. This can be calculated from the 
approach of Wylie and Sipe [511148] . The rate for a transition r — >■ s is given by 



kl 

where the Bohr frequency is Jujjrs = Er — Eg, and the thermal occupation number 
n{!jjrs) is evaluated at temperature T. The matrix elements of the dipole operator 
p, are written in Cartesian components, and G^.; is the electromagnetic Green 
tensor at the position z of the atom. Note that we normalize it such that eoGki 
has units of inverse volume. This formula also describes the absorption rate of 
thermal photons (uurs = —oJsr < 0) because 

[1 + n{-'UJsr)] 3m [Gki {z; -i^sr)] = n{ujsr) Jm [Gki {z; ujsr)] ■ (17) 

We get the total decay rate 7r by summing Eq. (116(1 over the final states js). Fig. [7] 
can be taken as an illustration of the terms in this sum because the electric field 
noise spectrum is proportional to 3m [Gj.; (2; ujrs)] [1 + n{ujrs)\ - The Green tensor is 
calculated for a gold surface using the formulas of Ref. [51 . The dependence on 
distance is shown in Fig.|9l Compared to free space (Planck spectrum, dashed line), 
transitions among Rydberg levels are significantly enhanced at short distances, 
while the rates oscillate in the opposite limit (distance comparable to the transition 
wavelength) due to the interference pattern mentioned above. Note the quite strong 
destructive interference for a transition dipole parallel to the surface (solid curve) 
which can be understood from the image dipole at a perfectly conducting surface. 

The transition rates scale with the electric dipole matrix elements of the Ryd- 
berg levels. We note that for a Bohr frequency h\ujrs\ >• An, the matrix elements 
are much smaller because a kind of radial selection rule suppresses changes in the 
principal quantum number by more than a few units (see also Refs. [361148] ). For 
a typical final state among adjacent Rydberg levels, for example An* < 2, we find 
an enhancement of the transition rate by a factor 0(5) at 2 = 10 fim compared 
to free space. The corresponding lifetime is reduced from ~ 150 //s (free space) to 
~ 30 /iS, estimating the matrix element by eagn*'^. We note that the rate for this 
generic pair of levels essentially exhausts the thermal decay rate summed over all 
final states, including photoionization [30]. We therefore expect a total lifetime 
somewhat below the ^ 30 /is figure estimated above from a single bound-bound 
transition. While this is still three orders of magnitude longer than the Rydberg 
excitation pulse, it may pose a serious challenge to gate errors below the fault tol- 
erance threshold of 10^*^. As a preliminary check, we have performed (see Table[2]) 
a calculation of the phase gate error when a finite lifetime 20 /is of the Rydberg 
level is included: an increase of the order of 10^'^ is indeed found. 

The radiation field also induces a van der Waals-Casimir-Polder shift on the 
Rydberg levels. We can estimate this in the London limit (transition wavelength 
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large compared to the atom-surface distance) since the Rydberg spacing A^g cor- 
responds to wavelengths in the centimeter range, 

^^'^^(^^ = 8^ife^ • 

The expectation value of the squared dipole is of the order of ^{eagii*'^)'^ (see, 
e.g., Ref. [54j l and gives a shift ^ 1.6 MHz at 10/im, consistent with the findings 
of Ref. [3S] where the electric quadrupole contribution is analyzed as well. This is 
not far from the dipole-dipole interaction (as it must from the scaling), but still 
small enough not to perturb it. The level shift changes only weakly across the 
levels shown in Fig. |6] and does not induce significant detunings. Its main impact 
is therefore to pull the Rydberg atoms towards the chip during the pulse. This 
effect which excites motional states in the trap, could be compensated for by the 
exciting laser pulse in a similar way as the momentum exchange between the two 
Rydberg atoms (see Sec. I2.3|l . 

We finally turn to the question how thermal radiation is shifting the Rydberg 
levels. This could be significant since even room-temperature blackbody radia- 
tion produces a sizable electric field above 100 V/m, albeit over a wide frequency 
range. It has been pointed out that the Casimir-Polder potential is essentially 
temperature-independent, due to cancellations between different transitions, on 
the one hand, and virtual and real photon exchange, on the other [l9]. This holds 
provided the typical transition wavelengths are large compared to the atom-surface 
distance, which is indeed the case for a Rydberg atom. We have checked that this 
result can be understood in a simple way starting from the dynamical polariz- 
ability of a free electron and integrating over the thermal radiation spectrum. For 
an analysis at zero temperature, see Ref. [55 and references therein. The Planck 
spectrum gives a free-space level shift, common to all weakly bound Rydberg lev- 
els, of the order of afs{kBT)^ / {nieC^) = 0(2 kHz), where offg is the fine structure 
constant and mgC^ the electron's rest energy. This is consistent with the value 
quoted in Ref. [36_. Near the chip surface, the shift is modified by a factor \t/z, 
where Ay = hc/ksT « 7.6 /im is Wien's thermal wavelength. The result is still 
negligible on the energy scale set by the Rydberg dipole-dipole interaction, and 
we provide a detailed discussion elsewhere. 



4 Conclusions 

We have used optimal control theory to determine the shortest gate duration for 
a controlled phase gate based on resonant excitation of neutral rubidium atoms to 
Rydberg levels that show a long-range dipole-dipole interaction. The parameters 
were chosen to be similar to those of the experiment reported in Ref. [26] , in partic- 
ular a distance between the atoms of 4 ^m, a trapping frequency of aj/27r ^ 276 kHz, 
and near-resonant two-photon interaction to the Rydberg level via the 5pi/2 level. 
It turns out that in this setting the fidelity of the controlled phase gate is limited to 
a few percent by the strong decay from the intermediate level, despite suppressing 
the population of the intermediate level as much as possible. The optimal pulses 
correspond to STIRAP-like transitions to and from the Rydberg state, i.e. an adia- 
batic solution to the control problem is found. Neglecting spontaneous decay from 
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the intermediate level (which would correspond to a different excitation scheme to 
the Rydberg level), gate errors of the order of 10~^ are obtained for T = 30 ns and 
larger. This lower limit on duration is essentially due to the interaction strength 
of the two atoms in the Rydberg state. The control pulses induce a sequence of 
Rabi flops, and the ensuing nuclear dynamics is strongly non-adiabatic. Further 
reduction of the gate error to values below the fault tolerance threshold of lO"'' 
requires the calculation to include the trapping potential for the qubit states dur- 
ing the gate. Otherwise, the motional state of the atoms at the end of the gate is 
not well-defined to this precision, and higher lying trap states get excited. 

These conclusions hold in a general setting where our parameters for trap sepa- 
ration and frequency are applicable. In order to estimate whether a Rydberg phase 
gate can be implemented on an atom chip to yield a universal quantum computer 
in a scalable setting, all relevant noise sources specific to the chip environment 
need to be considered. To summarize the estimates on sensitivity to noise due to 
the chip environment, we have found serious issues for Rydberg atoms held at 
a distance of the order of 10 /im, due to linear and quadratic Stark shifts. The 
reason is the contamination by impurity atoms of the chip surface at densities 
higher than 100 ^m"'^, still much less than a monolayer. This would create, above 
a metallic surface, electric fields exceeding 0(1 V/m), and reduce significantly the 
strong dipole-dipole interaction between the Rydberg atoms. This slows down the 
quantum gate, and a sizable dephasing rate arises from the beating between stray 
fields and fiuctuating patch potentials, in particular for the highly polarizable 56f 
state. Operating the chip at lower temperatures would reduce the patch charge 
noise |431I45| . The radiative decay of the Rydberg state is also enhanced by nearly 
an order of magnitude compared to free space at zero temperature: we estimate 
a lifetime in the 20 ... 30 /^s range. As it stands, this could provide a fundamen- 
tal lower limit around 10"'^ to the gate error. Possible improvements may exploit 
a Purcell effect to suppress, by destructive interference, radiative transitions for 
certain orientations of the transition dipole, similar to the suggestion in Ref. [53J. 
This may be achieved with suitably polarized laser pulses. We therefore estimate 
that reaching errors below the fault tolerance threshold for a Rydberg phase gate 
on an atom chip is challenging but possible. 
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